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ABSTRACT 


The  effect  of  the  position  of  the  X-ray  unresolved  K 


- >^a2 


doublets,  diffracted  from  powder  specimens,  on  the  precision  of  the 
calculated  lattice  constant  has  been  determined  using  a least-squares 
analysis.  . An  analytical  procedure  to  synthesize  CuK^  doublet  X-ray  dif- 
fraction peaks  with  X-ray  characteristic  lines  (half  widths  ranging  from 
0.1°  to  0.4°  29)  has  been  applied  in  correcting  the  weighted  wavelength 
of  the  doublet  peak  position.  A series  of  correction  curves  was  estab- 
lished from  which  the  true  26  peak  position  of  the  weighted  wavelength 
could  be  determined  from  the  measured  29  peak  position. 

Data  taken  from  silicon  powder  and  an  yttrium  aluminum  garnet  (YAG) 
powder  standard  to  calculate  the  lattice  constants  indicate: 

a.  A loss  in  precision  of  approximately  one  part  in  25,000  in 
determining  the  lattice  constant  is  obtained  when  calculated  using  dou- 
blet 26  peak  positions  applied  in  the  conventional  manner  (X  weighted  = 
1.54178  A,  at  peak,  for  copper). 


b.  An  increase  in  precision  is  obtained  by  using  low  as  well  as 
high  26  angle  diffraction  peaks  for  calculating  the  lattice  constant. 

The  effect  of  the  measured  characteristic  line  width,  simulated  by 
modifying  the  slit  system,  on  the  calculated  lattice  constant  was  deter- 
mined in  the  case  of  the  YAG  powder  specimen. 

\ 

Suggestions  to  increase  the  precision  of  the  lattice  constants  de- 
termined from  X-ray  diffraction  powder  samples  are  presented. 
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The  determination  of  the  lattice  constant  ao  of  a unit  cell,  using 
the  positions  of  characteristic  X-ray  lines  diffracted  from  a powder  sam- 
ple, has  been  treated  by  many  investigators.  Calculations  of  the  lattice 
constant  made  by  extrapolating  d values  against  various  functions  of  the 
Bragg  angle  0 have  been  made  using  a least-squares  analysis.  A summary 
of  the  systematic  error  functions  associated  with  an  X-ray  diffractometer 
has  been  presented  by  Vassamillet  and  King  (1),  who  found  that  the  best 
precision  to  be  expected  using  a commercial  diffractometer  is  approxi- 
mately 1:20,000. 

A comprehensive  procedure  that  has  been  developed  by  Mueller,  Heaton 
and  Miller  (2)  (MHM  method)  involves  varying  the  systematic  error  func- 
tions and  computing  the  change  in  d for  each  systematic  error  function, 
or  combination  of  functions,  until  the  error  in  the  computed  lattice 
constant  has  been  reduced  to  an  acceptable  level. 

In  a study  (3)  of  the  systematic  errors  involved  in  the  computation 
of  ag  for  yttrium  aluminum  garnet  powder  Y3AI5O12  (YAG) , it  was  found 
that  the  precision  of  the  computed  lattice  constant  was  a function  of 
the  following. 

(1)  The  slope  of  the  least-squares  straight  line  of  the  a^  values, 
where  a^  is  the  lattice  constant  computed  for  each  diffraction  line,  ver- 
sus cotSj,.  A slope  value  of  0 indicated  that  the  systematic  errors  had 
been  adequately  corrected  for  and  yielded  the  most  precise  ag  values. 

(2)  A statistical  error  term,  i.e.,  standard  error,  absolute  aver- 
age deviation,  rms  deviation  of  functions  of  aj^  and  sinSj^  from  their 
computed  functions  or  other  error  terms.  A comparison  of  these  error 
functions  ultimately  yield  the  same  lattice  constant.  IVhereas  a differ- 
ence was  seen  between  error  functions,  their  relative  values  were  useful 
in  establishing  the  same  most  precise  ag  values. 

(3)  The  number  and  selection  of  the  diffraction  peaks.  Using  all 
35  diffraction  peaks  (unresolved  K - K doublets  and  high-angle  re- 

solved  peaks)  yielded  the  highest  error.  Restricting  the  analysis  to 
the  12  highest  26  angle  peaks  resulted  in  a reduction  of  8%  in  error. 
Taking  the  20  highest  26  angle  peaks  increased  the  precision,  reducing 
the  error  28%.  The  fourth  case  gave  the  highest  precision  (42%  reduc- 
tion in  error  term),  and  occurred  when  23  high-  and  low-26  angle  peaks 
were  chosen,  eliminating  those  doublet  peaks  in  mid-26  range,  along  with 
the  diffracted  peaks  of  low  intensity. 


1 . VASSAMILLKT,  L.  I',  and  KING,  H.  W.  Precision  X-Ray  Diffractometry  Using  Powder  Specimens  in  Advances  in 
X-Ray  Analysis,  W.  M.  Mueller  el  al.,  ed„  v.  6,  Plenum  ftess.  New  York,  1963,  p.  142-157. 

2.  MUELLER.  M.  IL,  HEATON,  L„  and  MILLER,  K.  T.  Determination  of  Lattice  Parameters  with  the  Aid  of  a 
Computer.  Acta  Cryst.,  v,  13,  I960,  p.  828-829. 

3.  GAZZARA,  C.  P.  The  Kffect  of  Systematic  Errors  on  the  Measurement  of  Lattice  Constants.  Army  Materials 
and  Mechanics  Research  Center,  AMMRC  TR  70-29,  September  1970. 
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Recently,  in  an  attempt  to  relate  the  peak  height  to  the  integrated 
intensity  of  X-ray  diffracted  doublets  by  convoluting  the  characteris- 
tic peaks,  it  was  found  that  the  weighted  wavelength  approximation  was 
inaccurate  except  at  very  low  Bragg  angles  (4).  It  therefore  seemed  ap- 
propriate to  determine  what  the  effect  of  this  doublet  peak  position 
error  would  be  in  computing  the  lattice  constant  of  a material  whose  lat- 
tice constant  was  knowTi  and  to  afford  a means  for  correcting  for  the  dou- 
blet position. 


PROCEDURE 


A doublet  peak  position  correction  curve  could  be  generated  for  CuK,^ 
radiation  from  the  curve  shown  in  Figure  1 (4),  relating  the  value  of  the 
positional  parameter  Xj  of  the  analytical  expressions  describing  the  K 

^ 1 


and  lines  to  a,  the  function  defined  as  the  peak  separation 

in  terms  of  the  wavelength  dispersion.  Note  that  it  is  only  in  the  lin- 
ear portion  of  the  curve,  that  the  weighted  peak  position  approxima- 
tion is  valid  (i.e.,  Xj  = a/3). 


Figure  1.  Plot  of  posi- 
tional parameter  Xi 
versus  CuKq^  CuKu2 
diffracted  peak  separa- 
tion a,  assuming  Gaus- 
sian and  Gaussian  plus 
Cauchy  distribution. 


4.  CiAZZARA.  P.  Peak  Height  Approximation  for  X-Ray  Diffracted  Integrated  Intensity  in  Advances  in  X-Rjy 

Analysis,  R.  W.  (iould.  S.  Barrett.  J.  B Newkirk,  and  C.  O.  Ruud.  ed..  v.  19,  Kendall  Hunt,  Iowa,  1976,  p.  735-748. 
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The  solid  curve  in  F'igure  1 was  generated  assuming  a Gaussian  dis- 
tribution of  the  K„,  and  diffraction  peaks  as  a function  of  6.  The 

U \ Qt2 

dashed  curve  was  computed  by  fitting  a computed  characteristic  K peak 

to  an  experimental  CuK  (531)  Si  diffraction  peak,  where  the  upper  half 

was  assumed  to  be  Gaussian  and  the  lower  half  to  have  a Cauchy  distribu- 
tion (Gaussian  and  Cauchy). 

It  can  be  shown  that  the  equation  relating  the  change  in  wavelength 

to  a or  Xia  is:  ^ „ 

^ AX  = Xacjj^^/4.72  tan9 

where  is  the  characteristic  peak  width  at  half  peak  height.  Values 
of  can  be  obtained  from  Figure  2 from  the  measured  peak  width  at  half 

height  of  the  unresolved  peak  doublet  wy  . A series  of  curves  was 

-'2 

generated,  see  Figure  3,  relating  the  26  doublet  peak  position  to  the 
t correction  term  whereby 


26 


corrected 


= 26 


measured 


+ t, 


Figure  2.  ujhy^  versus  u>jy  curves  calculated  for 
various  levels  of  20  from  20°  to  55°  assuming 
Gaussian  and  Cauchy  peak  distributions. 


?9  degrees 


Figure  3 is  also  made  up  of  solid  lines  relating  to  Gaussian  dif- 
fraction peaks  while  the  dashed  curves  were  calculated  assuming  a Gaus- 
sian and  Cauchy  characteristic  diffraction  line. 

To  determine  the  effect  such  a t correction  would  have  on  the  lat- 
tice constant  precision,  CuK^^  diffraction  peak  positions  of  silicon  pow- 
der were  measured.  A CuK^^  powder  pattern  of  silicon  powder  with  a 
particle  site  of  less  than  30  microns  was  taken  with  a 0.05°  receiving 
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Figure  3.  Plots  of  peak  correction 
term  t versus  observed  diffracted 
20  peak  position  of  CuK^  dou- 
blets for  constant  values  of  why 


slit.  A scanning  speed  of  0.10°  (2e)/min  and  a 1°  (26)  beam  slit  was 
used  on  a G.E.  XRD-5  X-ray  diffractometer.  For  the  first  three  diffrac- 
tion peaks  an  value  of  approximately  0.13°  26  was  determined  from 

Figure  2 and  appfied  in  deriving  the  t correction  in  Figure  3.  A least"- 
squares  extrapolation  procedure  utilizing  a systematic  error  term,  cote, 
and  a weighting  factor,  tan6,  was  employed. 

The  plane  of  the  Si  powder,  characterized  as  a thin  specimen  (less 
than  0.002"  thick),  was  adjusted  with  respect  to  the  distance  perpendic- 
ular to  the  plane  of  the  specimen.  The  purpose  for  this  operation  was 
to  minimize  the  slope  (slope  = Aa^/AcotO^),  and  thereby  increase  the  pre- 
cision of  the  ag  measurement.  Values  of  a^^  are  shown  plotted  versus  6 
in  Figure  4 for  three  specimen  settings.  For  position  A,  which  repre- 
sents the  ideal  alignment  for  an  infinitely  thick  specimen,  the  extrap- 
olated value  of  ag  at  9 = 90°  is  too  high.  The  plotted  values  and 
computed  curves  applying  the  t correction  are  shown  below  the  uncorrected 
curves.  The  specimen  was  moved  inward  0.005"  to  position  B,  and  the  ef- 
fects of  the  least-squares  computation  are  shown  in  Figure  4.  Next,  the 
calculated  slope  values  for  positions  A and  B were  plotted  versus  ag  (see 
Figure  5)  and,  following  a known  procedure  (3),  the  correct  ag  value  was 
determined  for  slope  = 0.  Using  this  value  of  ag , position  C was  found 
by  varying  the  specimen  position  until  the  (111)  CuK^  peak  was  maximized 
at  the  proper  29  value  for  slope  = 0.  Plotting  values  of  a^  for  position 
C are  shown  in  Figure  4 and  the  respective  calculated  values  of  ag  are 
shown  in  Figure  5.  The  computed  results  of  this  procedure  are  summarized 
in  Table  1. 

Since  a prior  study  of  lattice  constant  measurements  was  performed 
with  YAG  powder,  it  was  next  deemed  appropriate  to  ascertain  the  effects 
of  the  t correction  on  a YAG  powder,  where: 


Table  1.  SILICON  LATTICE  CONSTANT  VALUES  MEASURED  FROM 
THIN  SPECIMEN  AS  A FUNCTION  OF  SPECIMEN  DISPLACEMENT 


(Powder  Size  Less  Than  30  ’pM,  0.05  Receiving  Slit, 
1"  Beam  Slit,  O.lO^/Min  Scan  Speed) 


Specimen 

Position 

Number  of 
Di ffracted 
Peaks 

t Correction 

ao,  A 

Slope 

Weighted 

RMS 

Deviation 

A (aligned 

12 

None 

5.4323 

-0.0038 

0.00130 

for  thick 

12 

Gaussian 

5.4325 

-0.0044 

.00120 

specimen) 

12 

Gaussian 

5.4325 

-0.0045 

.00124 

& Cauchy 

9 High  Angle 

- 

5.4322 

-0.0035 

.00083 

B (specimen 

12 

5.4310 

0.0043 

0.00083 

displaced 

12 

5.4312 

.0036 

.00070 

0.005" 

12 

5.4312 

.0036 

.00073 

i nwards ) 

9 High  Angle 

5.4310 

.0043 

.00028 

C (specimen 

12 

5.4315 

0.0008 

0.00072 

in  optimum 

12 

5.4317 

.0002 

.00049 

al i gnment) 

12 

5.4317 

.0002 

.00051 

9 High  Angle 

5.4316 

.0007 

.00034 

a.  many  more  diffraction  peaks  are  available  (12  Si  -+■  35  YAG)  ; 

b.  an  accurate  value  of  the  lattice  constant  was  available  for 
the  YAG  powder; 

c.  the  effect  of  uy,  and  the  lattice  constant  could  be  found  by 
varying  the  receiving  slit  width,  from  0.05°  26  to  0.10°  26  to  0.20°  26; 

d.  the  case  of  using  only  the  high  angle  diffraction  peaks  could 
be  examined;  and 

e.  unresolved  doublets  beyond  the  26  correction  range  could  be 
eliminated  (this  case  denoted  as  m) . For  a 0,05°  receiving  slit,  peak 
Zh^  = 52  was  eliminated;  for  a 0.10°  slit,  peaks  Zh^  = 52,  54,  56  were 
removed;  and  for  a 0.20°  slit  peaks  Zh^  = 52  to  88  were  taken  out  of 
the  least-squares  analysis.  The  results  of  the  analysis  of  the  thin 
YAG  powder  specimen  are  summarized  in  Table  2. 


RESULTS 

The  effect  of  the  error  in  assigning  a CuK,^  wavelength  of  1.54178  X 
to  the  unresolved  characteristic  line  doublet  peak  position,  on  the  cal- 
culated lattice  constant  of  either  silicon  or  YAG  always  lowers  the  true 
value  of  ag  by  approximately  one  part  in  25,000  (see  Figure  5).  It  can 
also  be  seen  from  Table  1 and  Figure  4 that  the  t correction  is  practi- 
cally the  same  whether  the  and  K„,  lines  are  considered  to  be  Gaus- 

cxj  oi2 

sian  or  Gaussian  and  Cauchy.  It  is  also  evident  that  with  a t correction, 
the  error  terms  decrease  and  the  slope  decreases  when  the  systematic  er- 
ror function  is  properly  treated.  Table  1 and  Figure  5 show  the  effect 
on  the  analysis  resulting  from  eliminating  the  first  three  CuK^^^  silicon 
powder  diffraction  peaks  and  using  the  remaining  nine  higher  angle  peaks. 


Position  A - Optimum  alignment  of  thick 
specimen 

Position  B • Sample  displaced  inward 
0.005  ' from  Position  A 
Position  C - Sample  aligned  maximizing 
diffraction  peak  Ih^  ■ 3 for  Oq 
determined  from  Figure  5,  for  slope  ■ I 


O CuKg  Diffraction  Peak  Doublets  Uncorrected 

□ CuKg  Diffraction  Peak  Doublets  Corrected 
for  Peak  Position  from  Figure  3 
(Gaussian) 

O CuKg  Diffraction  Peak  Doublets  Corrected 
for  Peak  Position  from  Figure  3 
(Gaussian  and  Cauchy) 


Error  Bar  Height  for  Precision 
of  1 0.005°  ( B ) 


Figure  4.  Values  of  lattice  constant  a;  calculated  from  diffracted  (oeak  position  of  family 
of  planes  {hkl},  versus  Bragg  angle  0 for  a thin  silicon  powder  specimen. 


Q With  CuKg  Doublet  Correction 
O Without  CuK 

A \ine  High  Angle  ResolvedCuK^ 
K Diffraction  Peaks 


Doublet  Correction 


Figure  5.  Lattice  con- 
stant ao  versus  slope 
term  for  thin  silicon 
powder  specimen  show 
ing  effect  of  t correc- 
tion and  the  use  of 
high  angle  diffraction 
peaks  only. 
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The  data  points  in  Figure  5 are  approximately  the  same  as  those  obtained 
using  all  the  diffraction  peaks,  without  the  t correction,  in  agreement 
with  earlier  work  (3). 

From  Figure  4,  notice  the  effect  of  improperly  treating  the  system- 
atic error  in  extrapolating  to  ag.  In  this  case  an  error  of  one  part  in 
5000  in  ao  is  possible,  indicative  of  the  importance  of  considering  the 
slope  term,  unless  a procedure  such  as  the  MIM  method  is  employed. 

The  lattice  constants  computed  from  the  experimental  data  obtained 
from  the  YAG  powder  may  be  seen  in  Table  2 and  in  the  plot  of  Figure  6. 
The  value  of  the  lattice  constant  at  slope  = 0 was  determined  from  an 
earlier  NWM  analysis  to  be  12.0107  X (3). 

In  the  case  of  YAG,  it  can  be  seen  that  the  most  precise  values  of 
ao  are  obtained  with  a 0.05°  receiving  slit,  progressing  from  using  the 
uncorrected  doublet  peaks,  to  employing  the  t correction,  to  using 
the  high  angle  peaks  alone.  These  results  represent  the  only  conditions 
where  the  high  angle  data  gave  the  best  results,  in  conflict  with  ear- 
lier work  (3) . The  lowest  lattice  constant  values  are  obtained  with  the 
intermediate  size  receiving  slit  (0.10°)  when  only  the  high  angle  peaks 
are  utilized.  For  both  cases,  however,  the  change  in  lattice  constant 
with  slope  moves  along  the  same  line  shown  in  Figure  6.  For  the  0.10° 
receiving  slit,  when  eight  additional  mid-29  angle  diffraction  peaks  are 
removed  from  the  analysis,  the  change  in  lattice  constant  with  slope  de- 
parts from  the  curve.  For  the  0.20°  receiving  slit,  the  slope  values 
increase,  indicating  a further  change  in  the  magnitude  of  the  systematic 


Table  2.  YAG  LATTICE  CONSTAfJT  VALUES  MEASURED 
AS  A FUNCTION  OF  RECEIVING  SLIT  SIZE 


(Powder  Size  Less  Than  30  „M,  1"  Beam  Slit, 
0.10'’/Min  Scan  Speed) 


Recei v- 
ing 
SI  i t 
(2"°) 

Number 

of 

Di ffracted 
Peaks 

t Correction 

O 

, A 

SI  ope 

Weighted 

RMS 

Devi  ation 

0.05 

35 

None 

12.0068 

0.0063 

0.00347 

34(m) 

None 

12.0069 

.0060 

.00330 

35 

Gaussian 

12.0071 

.0056 

.00329 

34(m) 

Gaussian 

12.0072 

.0053 

.00305 

23  High 

12.0075 

.0047 

.00195 

Angle 

0.10 

35 

12.0052 

0.0081 

0.00380 

32(m) 

12.0054 

.0076 

.00357 

35 

12.0055 

.0075 

.00401 

32(m) 

12.0057 

.0070 

.00371 

21  High 

12.0050 

.0085 

.00233 

Angle 

0.20 

35 

12.0064 

0.0126 

0.00417 

26(m) 

12.0068 

.0113 

.00277 

35 

12.0066 

.0123 

.00428 

26(m) 

12.0070 

.0108 

.00280 

15  High 

12.0068 

.0113 

.00187 

Angle 
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errors,  with  a displacement  of  the  points  from  the  curve.  The  error 
terms  shown  in  Table  2 indicate  no  unusual  behavior  between  cases  except 
that  in  treating  the  0.20°  receiving  slit  case,  the  sensitivity  of  remov- 
ing the  mid-range  diffraction  peaks  seems  to  be  more  severe.  This  is 
expected  when  one  considers  that  the  doublet  correction  t has  not  been 
taken  into  account,  siiice  these  corrections  occur  beyond  the  angular 
range  considered  in  Figure  5. 


CONCLUSIONS 

1.  \ correction  for  the  26  X-ray  diffracted  peak  position  of  CuK^j 
unresolved  doublets  results  in  an  increase  in  lattice  constants,  calcu- 
lated from  a least-squares  extrapolation  procedure,  of  one  part  in  25,000 
Tor  silicon  and  YAG  powder. 

2.  The  same  peak  position  of  the  unresolved  doublet  results  from 

assuming  either  a Gaussian  or  a Gaussian  plus  Cauchy  distribution  of  the 
CuKj^  and  CuK^^j  X-ray  diffraction  lines. 


Figure  6.  Lattice  constant  ap  versus  slope  term  for  thin  YAG  powder  specimen 
showing  the  effect  of  t correction  and  tfie  use  of  high  angle  diffraction  peaks  only. 


3.  In  the  absence  of  an  extensive  procedure  such  as  that  of  MUM  (2) , 
an  extrapolation  procedure  utilizing  a plot  of  ag  versus  slope  can  be 
applied  to  minimize  the  effect  of  systematic  errors  and  increase  the 
precision  in  ag. 
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4.  The  most  precise  values  of  ao  are  usually  obtained  when  high- 
and  low-20  angle  diffraction  peaks  are  employed  in  a least-squares  analy- 
sis. This  is  particularly  true  when  few  diffraction  peaks  are  available, 
i.e.,  silicon  powder  with  CuK^j.  High  precision  can  also  be  achieved  when 
many  diffraction  peaks  are  available,  i.e.,  YAG  case,  utilizing  only  the 
high  20  angle  diffraction  peaks. 
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